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ABSTRACT: Hydrodynamic forces may significantly affect the motion of polymers. In sheet-like 
cavities, such as the cell's cytoplasm and micro-fluidic channels, the hydrodynamic forces are long- 
range. It is therefore expected that that hydrodynamic interactions will dominate the motion of polymers 
in sheets and will be manifested by Zimm-like scaling. Quite the opposite, we note here that although 
the hydrodynamic forces are long-range their overall effect on the motion of polymers vanishes due to 
the symmetry of the two-dimensional flow. As a result, the predicted scaling of experimental 
observables such as the diffusion coefficient or the rotational diffusion time is Rouse-like, in accord 
with recent experiments. The effective screening validates the use of the non-interacting blobs picture 
for polymers confined in a sheet. 

Introduction 

When a polymer moves around inside a solution it drags the surrounding solvent and thus may 
mediate hydrodynamic forces between its monomers 1 . These forces are known to dramatically change 
the dynamic behavior of polymers as manifested by experimental observables such as the diffusion 

2 3 

coefficient D, the rotational diffusion time x R or the dynamic scattering spectrum ' . Extensive studies 
of the hydrodynamic interactions have distinguished two asymptotic regimes: In the Rouse regime 4 , the 
movement of each monomer is determined solely by the local forces. In this regime the scaling of D and 
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t r as the function of the polymer's molecular weight M is D ~ M ' and z R ~ M .At the other limit, the 
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Zimm regime 5 , the many-body hydrodynamic interactions between the monomers dominate their 

- 1 /2 3 /2 

motion. The scaling in this regime changes to D ~ M~ and tr ~ M , and may differ from the Rouse 
scaling by several orders-of magnitude. Hence, beside the theoretical interest, the question of whether 
hydrodynamic forces affect polymer dynamics has a vast practical importance in a variety of 
experimental setups and engineered devices where polymers are transported and manipulated ' . 

The present note discusses the effect of the hydrodynamic forces on polymer solutions that are 
confined in a sheet, a geometry where two of the dimensions are much larger than the third, the height h. 
This configuration is sometimes termed slit, crack or 'quasi two-dimensional' 7 ' 8 . The effect of the 
boundaries is expected to be significant when the sheet's height is smaller than the gyration radius of the 
polymer. Such geometry is typical to biopolymers confined in living cells or to macromolecules flowing 
through the cracks of a porous media. Flat channels are also the basic wiring of the rapidly developing 
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micro-fabricated fluidic devices ' and their design therefore relies on the understanding of quasi-2D 
hydrodynamics. 

A series of recent experimental and theoretical studies have demonstrated that colloidal beads 
suspended in a fluid sheet between two plates exhibit anomalously long-range hydrodynamic 
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interactions " . Since the basic dynamic properties of polymers can be represented by a set of colloidal 
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beads connected along a chain ' , it is natural to expect that the same long-range forces also dominate 
the dynamics of a quasi-2D polymer solution and lead to Zimm-like scaling. In contrast, the simple 
conclusion of this note is that although the hydrodynamic forces are indeed long-range their overall 
effect vanishes by symmetry. As an immediate result, we predict Rouse-like scaling for polymers 
confined in flat channels or sheets, e.g. D ~ M' 1 and tr ~ M 2 , in accord with recent experiments 14 and 
simulation studies 15 ' 16 . 

In what follows, we first present a simplified derivation of the hydrodynamic forces in sheets. We then 
examine the impact of these forces on polymer dynamics and conclude that their net effect vanishes. 
Finally, we discuss the validity of this result to geometries where the boundaries are not parallel and its 
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relevance to the blob picture ' of 2D polymer dynamics. 
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Results and Discussion 

Quasi-2D hydrodynamics at low Reynolds number. Polymers flow in a sheet typically at a very low 
Reynolds number and the hydrodynamics is therefore that of the viscous Stokes flow. It follows from the 
linearity of the Stokes flow that the velocities u n of particles suspended in a viscous solution are linearly 
related to the forces acting on them f n by the Oseen tensor H 19 , 

(1) u„=XH(r„-r m ).f m , 

m 

where for the present 2D geometry H is a 2x2 tensor that depends on the distance between the particles 
J"ji _ r m . 

To examine the hydrodynamic forces in quasi-2D geometry we consider the flow-field around a coin- 
like disk of radius a that moves horizontally at a velocity U in the middle of a sheet bounded by the z = 
±h/2 planes (figure 1). In principle, given the no-slip boundary conditions it is straightforward though 
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rather cumbersome to derive the Oseen tensor from Stokes equations ' . Nonetheless, the equations 
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can be further simplified by employing the lubrication approximation used in the Hele-Shaw cell ' : 
Since the height h is much smaller than the horizontal dimensions, it restricts the profile of the velocity 
v to be approximately that of the parabolic Poiseuille flow, 



\(x,y,z) = u(x,y) 

(2) 
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u(x,y) = -—vp(x,y), 

where rj is the viscosity of the host fluid and the mid-plane velocity \\{x,y) is proportional to the gradient 
of the pressure p. Equation 2, termed sometimes the Darcy approximation, is equivalent to a 2D 
potential flow, u(x,y) = V§ with the effective potential § = -h p/8rj. Accounting for incompressibility, 
V-v = 0, one finds the Laplace equation, A(|) = 0. 
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The flow around the disk is found by solving the Laplace equation with the boundary condition of zero 
mass-flux through the edge of the disk. This yields a potential and a velocity field of a two-dimensional 
mass-dipole of strength proportional to the area and velocity of the disk, d = a U, 



where f is a unit vector in the direction of the position r. Summing the momentum flux through the 
disk's edge we find the overall drag exerted by the host fluid, which is counterbalanced by the external 
force f required to move the disk, f ~ r\d/h. And so, exerting a force f on a particle induces at r a velocity 
u given by (3), which in turn may drag other particles. This relation between the induced flow and the 
force determines the quasi-2D Oseen tensor (1), 



This dependence on distance and direction of the hydrodynamic interaction was derived in exact 
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solutions for point forces ('stokeslets') confined between parallel plates ' . Nevertheless, the 
simplified derivation adds some insight by clarifying the relation to 2D potential flows. The dipolar flow 
(3) approximates well the far-field velocity induced by the monomers of a polymer in a sheet. 

The effect of hydrodynamic forces on polymers in a sheet. To explore the basic dynamic properties 
of a polymer, its monomers are represented as a set of Brownian particles connected along a chain by 
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elastic springs . This classical picture is expressed by the generalized Langevin equation for the 
monomer positions r„ l , 
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(4) H(r)~-4-(2f®f-l). 
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Figure 1. (a) Illustration of the flow field around a thin disk of radius a moving in a uniform sheet at a velocity 
U. The streamlines and the potential § are that of a two-dimensional dipole of strength d = a U. It is evident that 
the orientational average of the velocity vanishes due to symmetry, (b) The same disk moving at the same 
velocity in a sheet of height that varies as h(r) = ho (1 + e-r), with a slope |£|=l/(5a). Although the streamlines are 
significantly distorted the orientational average of the perturbed velocity ( 9) still vanishes. 
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This may be derived from (1) by noting that the overall force acting on the m-th monomer is the sum of 
the elastic force, which is the gradient of the elastic energy E, and the stochastic thermal force f m . Zimm 
has devised a way to simplify the analysis of the nonlinear equation (5) by 'pre-averaging' the Oseen 
tensor (4) over the equilibrium distribution of distances between the monomers r nm = r n - r m 5 , 

(6) H(r„ m )^(H(r„J)D ; 2r nm ®r nm -l) , 

where we used the fact that the ensemble averages over orientations and distances are independent. 

The equilibrium average of the monomer-to-monomer separation scales as (\r nm \ 2 ) ~ \m-n\' 2v , with 
Flory's exponent v ~ 3/(D+2), where v = 3/4 in 2D and v = 3/5 in 3D. This algebraic decay of the 
average indicates strong hydrodynamic forces (Ref. 1 p. 93) and one might expect to observe the Zimm 
dynamics. However, due to the unusual symmetry of the Oseen tensor, its average over orientations 
vanishes, (2r nm ®r nm — l) = . This immediately implies that the relevant hydrodynamic forces acting on 
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a monomer are only due to the nearby monomers in the local 'blob' of radius h around the monomer ' 

24 

. All the other, long-range hydrodynamic interactions may be disregarded. It follows directly that 
quasi-2D polymer solutions are always in the Rouse regime. We note that this conclusion relies on the 
accuracy of the pre-averaging approximation . This result is in contrast to the situation in 3D polymer 
solutions where determination of the Rouse and Zimm regimes and the crossover between them is a 
matter of subtle modeling and observation . 

The physical grounds that make the long-range hydrodynamic forces insignificant for the motion of 

quasi-2D polymers can be understood by re-examining the flow field in (3). The average of the Oseen 

tensor over all orientations (6) is actually the average over the dipole field that must vanish due to the 

symmetry of the dipole. As pointed out by the authors of references 8-10 this is very different from the 

situation in 3D Stokes flows that exhibit a momentum monopole. There, the pre-averaged contribution 

of the momentum monopole does not vanish and leads to the familiar Zimm scaling. In quasi-2D flows 
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the monopole contribution decays exponentially because the confining plates absorb the momentum and 
the remaining leading term is the mass dipole. 

Screening of hydrodynamic forces in non-uniform sheets. To explore the extent of validity of our 
conclusion that the overall effect of the hydrodynamic forces vanishes due to symmetry, we examine a 
geometry where the height of the fluid sheet is non-uniform (figure lb). If the height h changes 
moderately around an average ho, we can linearly approximate it as h{r) ~ ho (1 + £-r), where the 
gradient is Vh ~ hot. This may describe a sharp wedge or a slit. The slope is moderate enough to keep 
the Darcy approximation (2) valid. However, since the velocity is proportional everywhere to h and the 
flux to h (integrating along the z-direction) the Laplace equation is modified into 



with the effective potential § = -ho p/8rj. Substituting the linear dependence of h, equation 7 can be 
rewritten as (1 + 3s-r)A(|) = - 3s -V<|> . 

Since the gradient £ is small we can expand the potential around the dipole flow of a uniform sheet § Q 
(3), § ~ §o + <h- The perturbation §\ may therefore be found from Acj>i = - 3s -V^o. Solution of the latter 
equation with the appropriate boundary condition on the disk gives 



(7) V-((h/h fv</>) 



= 0, 



(8) 




The corresponding perturbation ui of the velocity, u = {hi ho) V§ = V§o + Ui , is 



(9) u 1 = 



2(£-r)V4 + V$ =— (ll(£-f)(d-f)f-(£-f)d-3(d-r)£) + o(£ 2 ). 
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The resulting perturbation of the flow field is shown in figure lb. Despite the significant distortion of 
the field, it is straightforward to verify that the average over all orientations of (9) vanishes, (ui)=0. 
Thus, up to first order in the gradient s, the pre-averaged hydrodynamic interaction vanishes in 
moderately non-uniform fluid sheets. This implies that the hydrodynamic interactions of a polymer are 
effectively screened also in non-uniform sheets if the height h does not change much over the polymer's 
radius of gyration R & \e\R g « 1. 

Conclusions 

The conclusion that the hydrodynamic forces are screened due to symmetry relies only on the 
soundness of the Darcy approximation. It will therefore remain valid even if the sheet is curved, for 
example into a spherical or cylindrical shell, as long as the height does not vary too steeply. 

The dipolar flow described here is not unique to polymer solutions in sheets. It was observed in a 
number of two-dimensional configurations 26 ~ 28 . Screening of long-range hydrodynamic forces was 
noted in electrophoresis and for Brownian particles in a thin film , where the momentum flux is 
dissipated by a phenomenological friction force. It can be expected to observe similar screening effects 
in other configurations whenever a mechanism that dissipates momentum is present. 
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Finally, we remark that the present effective screening justifies the use of the blob picture to derive 
the scaling of the diffusion coefficient and other related observables of polymer solutions in thin sheets 
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' . This is because the blob scaling relies on the fact that blobs do not exert hydrodynamic forces on 
each other due to the strong screening. The total friction of the polymer is therefore proportional to the 
number of blobs Nb times the friction of a single blob, which by itself is proportional to the blob size h. 
Taking into account the dependence of Nb on h and on the polymer's mass M, the blob scaling of the 
diffusion coefficient is Rouse-like, D ~ (Nbh)' 1 ~ h lv ' ! M A . This scaling of the diffusion coefficient is 
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well-known for polymers confined in a tube , where the hydrodynamic forces decay exponentially and 
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the polymer therefore exhibits Rouse-like dynamics ' . The present study suggests that the same 

scaling holds also for polymers confined in a thin sheet, although for very different reasons. Unlike the 
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flow in the tube, the long-range hydrodynamic forces between monomers of a polymer confined in a 
sheet are not screened, yet their overall effect vanishes due to symmetry of the flow field. 

The blob scaling agrees with recent diffusivity measurements of double-stranded DNA molecules 
confined in a slit-like channel 14 , where the blob scaling is presented in the equivalent normalized form 
DIDbuik ~ (Rglh) 1 ' 1/v . Nevertheless, the current measurement has reached only moderate values of the 
confinement parameter, R g /h ~ 0.25-0.9. In the high confinement regime of thinner slits and longer 
polymers, the blob scaling is supported only by simulation studies 15 ' 16 and is yet to be experimentally 
tested. We expect the h 1/v ~ ! dependence of the blob scaling to break down for ultra-thin sheets 29 whose 
height is of the order of a few persistence lengths. The M" 1 dependence merely manifests the locality of 
the forces and is therefore expected to remain valid even for very thin sheets . 



9 



References 



1. Doi, M.; Edwards, S. F., The theory of polymer dynamics. Oxford: 1986. 

2. Dubois- Violette, E.; de Gennes, P. G. Physics 1967, (3), 181. 

3. Shusterman, R.; Alon, S.; Gavrinyov, T.; Krichevsky, O. Phys Rev Lett 2004, 92, (4), 048303. 

4. Rouse, P. E.JChem Phys 1953, 21, (7), 1272-1280. 

5. Zimm, B. H. J Chem Phys 1956, 24, (2), 269-278. 

6. Squires, T. M.; Quake, S. R. Rev Mod Phys 2005, 77, (3), 977-1026. 

7. Darnton, N.; Bakajin, O.; Huang, R.; North, B.; Tegenfeldt, J. O.; Cox, E. C; Sturm, J.; Austin, 
R. H. J Phys: Cond Mat 2001, 13, (21), 4891-4902. 

8. Diamant, H.; Cui, B.; Lin, B. J Phys: Cond Mat 2005, 17, (31), S2787-S2793. 

9. Cui, B. X.; Diamant, PL; Lin, B. PL; Rice, S. A. Phys Rev Lett 2004, 92, (25), 258301. 

10. Diamant, PL; Cui, B.; Lin, B.; Rice, S. A. J Phys: Cond Mat 2005, 17, (49), S4047-S4058 

1 1 . Falck, E.; Lahtinen, J. M.; Vattulainen, I.; Ala-Nissila, P. Eur Phys JE 2004, 13, (3), 267-275. 

12. Alvarez, A.; Soto, R. Phys Fluids 2005, 17, (9). 

13. Riseman, J.; Kirkwood, J. G.J Chem Phys 1950, 18, (4), 512-516. 

14. Chen, Y. L.; Graham, M. D.; Pablo, J. J. d.; Randall, G. C; Gupta, M.; Doyle, P. S. Phys Rev E 
2004, 70, (6), 060901. 

15. Milchev, A.; Binder, K. J de Physique II 1996, 6, (1), 21-31. 

16. Hagita, K.; Pakano, H. J de Physique IV 2000, 10, (P7), 305-308. 

17. Pincus, P. Macromolecules 1976, 9, (3), 386-388. 

18. Pincus, P. Macromolecules 1977, 10, (1), 210-213. 

19. Happel, J.; Brenner, PL, Low Reynolds number hydrodynamics. Kluwer: Boston, 1988; p 553. 

20. Liron, N.; Mochon, S. JEngMath 1976, 10, (4), 287-303. 

21. Staben, M. E.; Zinchenko, A. Z.; Davis, R. H. Phys Fluids 2003, 15, (6), 171 1-1733. 

22. Bensimon, D.; Kadanoff, L. P.; Liang, S. D.; Shraiman, B. I.; Tang, C. Rev Mod Phys 1986, 58, 
(4), 977-999. 

23. Daoud, M.; de Gennes, P. G. J de Physique 1977, 38, (1), 85-93. 

24. Brochard, F.; de Gennes, P. G. J Chem Phys 1977, 67, (1), 52-56. 

25. Wang, S. Q.; Freed, K. Y.J Chem Phys 1987, 86, (5), 3021-3031. 

26. Alder, B. L; Wainwright, T. E. Phys Rev A 1970, 1, (1), 18. 

27. Long, D.; Stone, H. A.; Ajdari, A. J Colloid Interface Sci 1999, 212, (2), 338-349. 

28. Seki, K.; Komura, S. Phys Rev E 1993, 47, (4), 2377-2383. 

29. Dufresne, E. R; Altman, D.; Grier, D. G. Europhys Lett 2001, 53, (2), 264-270. 

30. Maier, B.; Radler, J. O. Macromolecules 2000, 33, (19), 7185-7194. 



10 



